Nuclear enhancement of universal dynamics of high parton densities 
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We show that the enhancement of the saturation scale in large nuclei relative to the proton is 
significantly influenced by the effects of quantum evolution and the impact parameter dependence 
of dipole cross sections in high energy QCD. We demonstrate that there is a strong A dependence 
in diffractive deeply inelastic scatteringand discuss its sensitivity to the measurement of the recoil 
nucleus. 
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PACS numbers: 24.85.+p,13.60.Hb 

The properties of hadronic and nuclear wave functions 
at high energies are of great importance in understanding 
multi-particle production in QCD. Especially intriguing 
is the possibility that the small Feynman x components of 
these vifavefunctions demonstrate universal behavior that 
is insensitive to the details of hadron or nuclear structure 
in the (large x) fragmentation region. 

The specific nature of universal small x dynamics 
in QCD follows from the strong enhancement of gluon 
bremsstrahlung at small x leading to a rapid growth of 
the occupation number of a transverse momentum mode 
kj_ in the hadron or nuclear wavefunction. However, it 
can maximally be of order l/ag (where as is the QCD 
coupling constant) because of non-linear multi-parton ef- 
fects such as recombination and screening which deplete 
the gluon density at small x In particular, the oc- 
cupation number is maximal for modes with k± ^ Qs, 
where Qsix), appropriately called the saturation scale, is 
a scale generated by the multi-parton dynamics. For a 
probe with transverse resolution l/Q^, this scale is man- 
ifest in a universal scaling form of observables as a func- 
tion of Q/Qs in a wide kinematical range in x and Q^. 

In addition to the strong x dependence generated by 
gluon bremsstrahlung, the saturation scale Qs has a 
strong A dependence because of the Lorentz contraction, 
in the probe rest frame, of the nuclear parton density. For 
large enough A and small enough x, the saturation scale 
is larger than Aqcd, the fundamental soft scale of QCD. 
In this letter, we will discuss the A and x dependence 
of the saturation scale and some of its ramifications for 
hard diffraction in nuclei. 

A saturation scale arises naturally in the Color Glass 
Condensate (CGC) 0] description of universal proper- 
ties of hadron and nuclear wavefunctions at small x. The 
CGC, when apphed to Deeply Inelastic Scattering (DIS), 
results [3, Hi leading order in as, in the dipole pic- 
ture of DIS [5| , where the inclusive virtual photon hadron 
cross section is 
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virtual photon to produce a quark-anti-quark pair of 

size r = |rj^| and -^^^{r±,x,h±) denotes the dipole 
cross section for this pair to scatter off the target at 
an impact parameter b^. The former is well known 
from QED, while the latter represents the dynamics of 
QCD scattering at small x. A simple saturation model 
(known as the GBW model of the dipole cross sec- 
tion, parametrized as -grg^ = 2(1 — e~'^^'^^'V^^^^'^) where 
Qg p(a;) = (xq/x)^ GeV^, gives a good qualitative fit to 
the HERA inclusive cross section data for = 3 ■ 10^'' 
and A = 0.288. However, the model does not contain the 
bremsstrahlung limit of perturbative QCD (pQCD) that 
applies to small dipoles of size r <C 1/Qs(a;)- 

In the classical effective theory of the CGC, to lead- 
ing logarithmic accuracy, one can derive the dipole cross 
section [3| containing the right small r limit. This dipole 
cross section can be represented (see however [7|) 



dip 
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= 2 [1 - exp {-r^Fix, r)Tp{h^))] , (2) 



Here 'ill rp{r±, z,Q) represents the probability for 



where Tp(h±) is the impact parameter profile function 
in the proton, normalized as J (Ph±Tp{h±) = 1 and F 
is proportional to the DGLAP evolved gluon distribu- 
tion [8| 

Fix, r^) = TT^as (m^ + A/r^) xg (x, f^l + /{2N,). 

(3) 

The dipole cross section in Eq. Q was implemented in 
the impact parameter saturation model (IPsat) [9] where 
the parameters are fit to reproduce the HERA data on 
the inclusive structure function F2. 

In general, the dipole cross section can range from in 
the r — > color transparency limit to 2, the maximal uni- 
tarity bound. The saturation scale Qs characterizes the 
qualitative change between these regimes; we shall here 
define Qs as the solution of j^^'*' {x, = 1/Q^{x,h±)) = 

2(1-6^1/4) [2|. 

The IPsat dipole cross section in Eq. ^ is applica- 
ble when leading logarithms in Q^ dominate over lead- 
ing logarithms in a;. At very small x, quantum evolu- 
tion in the CGC [2] describing both the bremsstrahlung 
limit of linear small x evolution as well as nonlinear RG 
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Predictions for shadowing compared to NMC 



FIG. 1: Left: 

data. Center: predictions for 12_Ff "/118-F2^ compared to 
NMC data at a; = 0.0125. Right: Likewise for = 5 GeV^ 
as a function of x. 



evolution at high parton densities, combined with a re- 
ahstic 6-dependence, is better captured in the bCGC 
model 0, lill. Both the IPsat model and the bCGC 
model provide excellent fits to a wide range of HERA 
data for x < 0.01 
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12l |. We will now discuss the pos- 
sibility that DIS ofl^ nuclei can distinguish respectively 
between these "classical CGC" and "quantum CGC" mo- 
tivated models. 

A straightforward generalization of the dipole formal- 
ism to nuclei is to introduce the coordinates of the indi- 
vidual nucleons {bj,^}. One obtains in the IPsat model, 



d2b. 



= 2 



1 - e-'''-^(^'''^^^*=i'^f(''-L-''^') 



(4) 



where F is defined in Eq. ([3]). The positions of 
the nucleons {b^j} are distributed according to the 
Woods-Saxon distribution TA(bj_j). We denote the 
average of an observable O over {bj^,} by (0)^ = 

/nti d^huTAihu)0{{hu})- The average differen- 
tial dipole cross section is well approximated byQ 



d^b, 



1 - ^^^dip 



(5) 



where, for large A, the expression in parenthesis can be 
replaced by exp (^—^-^^^^^'^dip^ parameters of 

the model come from either fits of the model to ep-data 
or from the Woods-Saxon distributions; no additional pa- 
rameters are introduced for eA collisions. The same ex- 
ercise is repeated for the bCGC model. 

In Fig. [1] (left), we compare the prediction of the IP- 
sat and bCGC models with the experimental data [2^ 
on nuclear DIS from the NMC collaboration llJI)- Fig- 



ure [T] (right) shows that the x dependence of shadow- 
ing for fixed in the IPsat model is very fiat. This is 
because the best fit to ep-data in DGLAP-based dipole 
models [^, Q is given by a very weak x-dependence at 



the initial scale Hq. A stronger ^-dependence also for 
large dipoles, such as in the in the GBW or bCGC mod- 
els, gives a stronger x-dependence of shadowing at fixed 
Q^. As shown in Fig. [T] (center), both the IPsat and 
bCGC models predict strong Q^-dependence (at fixed x) 
for shadowing. It is this latter effect which is primar- 
ily responsible for the shadowing effect seen in the NMC 
data. Precision measurements of F2/AF2 would shed 
more light on the relative importance of and x evolu- 
tion in this regime. 

We now turn to a discussion of the A and x dependence 
of the saturation scale. In a simple GBW type model, 
inserting a 0-function impact parameter dependence into 



Eq. ([5]) yields the estimate ^ 



r>2 a2/3 



0.26^i/3q2^ foj. 27ri?2 „ 20 mb and Ra ~ I.IA^/^ fm. 
The smallness of Qsa/Qs,p' due to the constant factor 
~ 0.26 has sometimes been interpreted [H, [l^ as a 
weak nuclear enhancement of Qs- We will argue here 
that detailed considerations of QCD evolution and the 
6-dependence of the dipole cross section result in a sig- 
nificantly larger nuclear enhancement of Qs- 

The effect of QCD evolution on Qs,A in the IPsat nu- 
clear dipole cross section is from the DGLAP-like growth 
of the gluon distribution. The increase in the gluon den- 
sity with increasing and decreasing (dominant) dipole 
radius r causes Qs grow even faster as a function of A. 
This is seen qualitatively for two different nuclei, A and 
B (with A > B), in a, "smooth nucleus" approximation 
of Eq. (g]) whereby X]i=i ^p(b_L — b^J is replaced by 
ATA(b^). We obtain 



QIa _ A TA{h^) Fix, Q2 a) F{x, Q2 a) 



QIb 



BTB(bx) F{x,Ql 



(6) 



The sca ling violations in F imply that, as observed in 
Refs. [i,[i3, the growth of Qs is faster than A^^^. Also, 
because the increase of F with is faster for smaller x, 
the A-dependence of Qs is stronger for higher energies. In 
contrast, the dipole cross section in the bCGC model de- 
pends only on the "geometrical scaling" combination [2^ 
rQs{x) without DGLAP scaling violations and therefore 
does not have this particular nuclear enhancement 
Precise extraction of the A dependence of Qs will play an 
important role in distinguishing between "classical" and 
"quantum" evolution in the CGC. 

A careful evaluation shows that because the density 
profile in a nucleus is more uniform than that of the pro- 
ton, the saturation scales in nuclei decrease more slowly 
with b than in the proton. The dependence of the satu- 
ration scale on the impact parameter is plotted in Fig. [51 
The saturation scale in gold nuclei at the median impact 
parameter for the total cross section 6mod. is about 70% 
of the value at & = 0; in contrast, Qs^p(&mod.) is only 
^ 35% of the value at 6 = 0. 

The A dependence of the saturation scale for various 
X is shown in Fig. O for the IPsat model on the left and 
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FIG. 2; Impact parameter dependence of the saturation scale 
for p, Ca and Au at 2: = 0.001 and = 1 GeV^. Details in 
text. 
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FIG. 3: Saturation scale at 6 = (open symbols) and b — 
bmad. (filled symbols) as a function of A for different x. The 
saturation scale for the proton is shown at ^ = 1 and by the 
arrows on the right. Left: IPsat model. Right: bCGC model. 



the bCGC model on the right. Note that in the IPsat 
model, at small x, Qg A(^mcd.) for gold nuclei is nearly 
identical to A^^^ times the value for the proton. The cor- 
responding enhancement for 6 = is significantly smaller 
as anticipated. The nuclear enhancement in the bCGC 
model is nearly as large, showing that it owes, for the 
kinematic range studied, much more to the relative im- 
pact parameter profiles (see Fig. ^ than to differences in 
QCD evolution. Nevertheless, the stronger A dependence 
of Qs js^{bmcd.) in the IPsat model relative to the bCGC 
model, especially at the smallest x values, clearly illus- 
trates the differences in quantum evolution between the 
models. The factor of 200^/^ w 6 gives a huge "oomph" 
in the parton density of a nucleus relative to that of a 
proton; one requires a center of mass energy ~ 14 times 
larger in an e+p collider relative to an e+Au collider to 
obtain the same Qg A(^mcd.)(a;)■ 
We will now focus on some interesting qual itative fea- 
tures of hard diffraction off nuclei (see also (20|). For sim- 
plicity, we will consider only the IPsat model here. The 
contribution of qq dipoles [28| to the inclusive diffractive 
cross section can be expressed as 



dt 



1 

16^ 



d r 1 dz 



alp{x,r,A^), (7) 



Eq. ^ over t, reads [2 
<^E,T = \ d^rx dz 



d^b 



dcr, 



dip 



d2b 



(8) 



The diffractive slope at i = depends on the size of 
the system. For small t ^ —1/R\ one expects a very 
steep f-dependence ^ expjDii?^} (with D ^ 1). In our 
picture of the nucleus as a "lumpy" collection of partially 
overlapping nucleons (Eqs. (H]) and an interesting 

question is whether this lumpiness shows up as a proton- 
like tail ^ explD'f/ip)} of the t-distribution. 

If one requires that the nucleus stays completely intact, 
the average must be performed at the amplitude 
level, and da^ / dt falls off very rapidly as exp{Dii?^}. 
Measuring the intact recoil nucleus at such a small t ex- 
perimentally at a future electron ion collider [23j is chal- 
lenging. Considerable physical insight into the diffrac- 
tive process can be obtained in events where the nu- 
cleus breaks up into color neutral constituents without 
filling the rapidity gap between the qq dipole and the 
nuclear fragmentation region. Such events correspond to 
performing the average (•)f,^ over the cross section [2^ . 
Eq. ([8|), instead of the amplitude. The difference between 
the two averaging procedures can be significant with in- 
creasing values of i; the result for calcium nuclei is shown 
in Fig. m The t dependence of the proton is shown as 
well; as expected, the "break up" cross section for cal- 
cium converges to A times the proton cross section with 
increasing t. 

The difference between "no break up" and "break up" 
integrated cross sections can be seen in Fig. [5] where we 
plot, as a function of A for fixed x and Q^, the double 



ratio i?^ifj , defined as the ratio (crj; -I- a 



) 



from Eqs. ^ and fS]) for a nucleus divided by the same 
ratio for a proton. For light nuclei {A < 40), R^^g < 1 
before going well above unity for large A. This is because 
the diffractive qq cross section is dominated by smaller 
impact parameters than the inclusive cross section; at 
small impact parameters, the matter density in a proton 
is larger than in light nuclei. For large nuclei, especially 
in the "break up" case, the nuclear enhancement of the 
fraction of diffractive events can be quite significant, up 
to a 100% enhancement relative to the fraction for a pro- 
ton. 

We thank A. Caldwell, C. Marquet and L. Motyka for 
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U.S. Department of Energy. 



where t = — Aj^ and crdip(a:;, r, is the Fourier trans- 
form of the dipole cross section with respect to b_L . The 
total diffractive cross section, obtained by integrating 
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FIG. 5: Ratio of the fraction of qq diffractive events in rmclei 
to the fraction in a proton plotted versus atomic number A 
for fixed x and values. 
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Our definition is equivalent to the saturation scale in the 

GBW model for a Gaussian dipole cross section. Note the 

difference with the convention in Ref. 

Data available from the E665 collaboration differ from 

the NMC data for a similar kinematic range. 

See, however Ref. on "pre-asymptotic" violations of 

geometrical scaling. 

Another interesting possibility that running coupling 
log a;-evolution results in a depletion of the A-dependence 

of Qe [ii. 

The contribution of higher Fock states will be addressed 
in future work. 

A comparison of Eqs. ((TJ and ^ makes explicit the uni- 
tarity limit crlJ'y < aL,T/'2, which is saturated by a black 

disk, -^5^(bx) = 26{R—h), leading to the i-distribution 

~ 4j Ji(^^i?)|V(-ii?^) ~ l+tR^/A+0{t^). In contrast, 
an exponential t distribution ~ e^* (see e.g. [21|) corre- 
sponds to a Gaussian in b, where the unitarity limit, sat- 



urated by 



= 2e 



-h-' /{2B) 



is< 

T ^ (yL,T/^- Specify- 



ing the bx dependence of the dipole cross section simul- 
taneously fixes both the normalization go of the cross 
section and the diffractive slope B — din cr^/ dt|t=o. 
(oq/B — 8tt and ao/B = An for black disk and Gaus- 
sian respectively.) 



